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ABSTRACT 

The type lb supernova 2010O was recently discovered in the interacting starburst galaxy Arp 
299. We present an analysis of two archival Chandra X-ray observations of Arp 299, taken 
before the explosion and show that there is a transient X-ray source at a position consistent 
with the supernova. Due to the diffuse emission, the background is difficult to estimate. We 
estimate the flux of the transient from the difference of the two X-ray images and conclude 
that the transient can be described by a 0.225 keV black body with a luminosity of 2.5 ± 
0.7 x 10 39 erg s _1 for a distance of 41 Mpc. These properties put the transient in between the 
Galactic black hole binary XTE J 1550-564 and the ultra-luminous X-ray binaries NGC 1313 
X-l and X-2. The high level of X-ray variability associated with the active starburst makes it 
impossible to rule out a chance alignment. If the source is associated with the supernova, it 
suggests SN2010O is the explosion of the second star in a Wolf-Rayet X-ray binary, such as 
Cyg X-3, IC 10 X-l and NGC 300 X-l. 
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I INTRODUCTION 

Type lb supernovae lack hydrogen in their spectrum, but show 
helium lines, while Ic supernovae also lack helium. They are re- 
garded as the core-collapse explosions of massive stars that have 
lost their hydrogen envelopes and thus have become Wolf-Rayet 
stars ( [Gaskell et al.||1986) . These Wolf-Rayet stars can originate 
from the most massive stars that can remove their hydrogen enve- 
lope on and just after the main sequence in strong stellar winds, 
or from lower-mass stars via a binary interaction that removes the 
hydrogen envelope (e.g. Podsiadlowski et al. 1992, No moto et al.| 
[T^95l >. 

In the last decade a new way of linking supernovae with their 
progenitors has become available. The growing archive of high- 
resolution images has made it possible to detect the progenitors of 
type II (hydrogen rich) supernovae in pre-explosion optical images 
(see[Smartt 2009 for a review). No progenitors of type lb or Ic 
supernovae have been found, despite deep pre-explosion images 
for 10 of them. However, the relative frequency of Voile to type 

II supernovae of around 0.4 suggests that binary interactions play 
an important role, as for a standard initial mass function, there are 
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not enough very massive stars that could lose their envelope via a 
stellar wind ( [Smartt|2009] ). 

An alternative method for directly detecting supernova pro- 
genitors is to use archival images in other wave bands, such as X- 
ray data. We started a program to search for type la supernova pro- 
genitors in Chandra X-ray data and have found one likely progen- 
itor and derived five upper limits (Voss & Nelemans 2008 ; Roelofs 
|eTaL] 2008, Nelema n7eral.|2008) . For type la supernova progen- 
itors the argument to look for X-ray progenitors is the sugges- 
tion that supersoft X-ray sources may produce type la supernovae, 
based on the accreting white dwarf model (Whelan & Iben 1973, 
Nomoto 1982 ). For type lb and Ic supernovae, the binary progen- 
itor scenario suggests the possibility that the lb or Ic explosion is 
the second supernova in the binary and thus before the explosion 
may have been part of a high-mass X-ray binary. Indeed, the rela- 
tive frequency estimates of Podsiadlowski et al. ( 1992 their figure 
16) suggest that the likelihood for a lb to be the second explosion in 
the binary is at least as high as it being the first one. It is therefore 
useful to constrain the X-ray properties of the direct progenitors of 
lb and Ic supernovae. 

Even more tantalising evidence for such scenario comes from 
the known Wolf-Rayet star in the X-ray binary Cyg X-3 ( van Kerk- 
|wijk et al.|1996| ) and the recent discovery of two X-ray binaries in 
which a (massive) Wolf-Rayet star orbits a black hole, IC 10 X-l 
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and NGC 300 X- l (Clark & Crowther 2004 , B auer & Brandt|20041 
|Prestwich eraT]|2007[ |Carpano et aL]|2007| ). Given the short life 
times of massive Wolf-Rayet stars, a lb or Ic supernova explosion 
is expected within next few million years in these systems. 

In this paper we report the discovery of a transient soft X- 
ray source at the position of the recent lb supernovae 2010O that 
exploded in the interacting starburst galaxy Arp 299 (IC 694/NGC 
3690). In Sect. [2] we discuss supernova 2010O and its host galaxy, 
in Sect. [3] the X-ray, optical and radio observations that we used 
and in Sect. [4] the results of the Chandra analysis. In Sect. [5] we 
discuss the implications of the finding, the possibility of a chance 
alignment and the scope for future work. 



2 SUPERNOVA 2010O IN ARP 299 

SN20100 was discovered on Jan 24, 2010 in the course of the Puck- 
ett Observatory Supernova Search ( Newton & Puckett 2010) and a 
spectrum taken with the Nordic Optical telescope on Jan 28, 2010 
showed it to be a type lb supernova (M attila et al.|2010] >. The su- 
pernova resides in the Eastern part of the interacting galaxy Arp 
299, which is often referred to as IC 694 in the literature, desig- 
nating the Western part as NGC 3690. However, as this convention 
is somewhat confusing] we use the name Arp 299 to refer to the 
whole system and use Eastern and Western part to identify the two 
galaxies. 

Arp 299, at a distance of about 41 Mpc, is a pair of inter- 
acting galaxies giving rise to a spectacular starburst (e.g. [Gehrz 
et al. 1983) and has in the last 20 years produced seven supernovae: 
1992bu (type II), 1993G (type II), 1998T (type lb), 1999D (type 
II), 2005U (type II or lib), 2010O (type lb) and SN2010P (type un- 
known) or eight if including the possible radio supernova reported 
by |Huang et al.|(T990| . The star formation rate is estimated to be 
ranging from 22-140 M yr _1 |Alonso-Herrero et aL]|2000| ) and 
Chandra X-ray images reveal extended emission and at least 18 
discrete sources of which one is an obscured AGN in the Western 
part of the interacting galaxy (Zezas et al. 2003). Radio observa- 
tions show several discrete sources, as well as extended emission 
in the nuclei of both parts which have been resolved in 30 indi- 
vidual sources using VLBI (e.g.|Gil de Paz et al.|2007l|Neff et al.| 
|2004l |Ulvestad||2009| |Perez-Torres et al.||2009| >, likely all young 
supernova remnants. 

Clearly Arp 299 is a very complex stellar system with tremen- 
dous activity, making association of different sources with each 
other difficult. Nevertheless, we will analyse the pre-supernova X- 
ray data of Arp 299 to search for possible emission of the progenitor 
of2010O. 



3 OBSERVATIONS 

There are Chandra images at two epochs before the supernova 
(taken in 2001 and 2005) , 25 ks ACIS-I an d 10 ks ACIS-S. The 
2001 data is published in |Zezas et ah] ( [2003] ), while the 2005 DDT 
observation (aimed at detecting an X-ray flare from the AGN in the 



Table 1. List of pre-SN Chandra observations. The pointing offset between 
the source and the center of the Chandra field is given in the last column. 



1 |Hibbard & Yun|fl999j discuss that it is likely, although not conclusively 
so that Swift referred to the compact dwarf to the northwest of Arp 299 for 
IC 694. All supernovae in either Eastern or Western part of Arp 299 (in- 
cluding SN 2010O) are classified as "in NGC 3690" by the Central Bureau 
for Astronomical Telegrams of the IAU. 



Instrument ID 



date 



exp. time offset 



ACIS-I 
ACIS-S 



1641 

6227 



2001-07-13 
2005-02-14 



24850 
10320 



1.2' 
0.4' 
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Figure 1. Hubble Space Telescope ACS F814W image of the region around 
SN2010O. The three circles indicate (right to left): our best optical position, 
the reported discovery position (consistent with the position shown in Bond 
|et al.|(2010) ), and the best position of the X-ray transient (Sect. |4j. The 
radius of the circles is Of 5. 



Western part) has not been published. In table [T] we show a log of 
the observations. 

We checked the alignment of the two X-ray images by com- 
paring the positions of seven isolated sources in the two images 
resulting in a difference of Ax =0.02±0.05; Ay = -0.11 ±0.07 pix- 
els. With a pixel scale of 0"491 this translates into a shift of 0"05 
with marginal significance. We therefore chose not to correct for 
this possible shift. 

We checked the absolute astrometric accuracy of the X-ray im- 
ages by comparing them to USNO B1.0 and 2MASS positions, but 
there are only two matches and they are off- axis in the X-ray im- 
ages so we can only conclude that they agree to within 1 '.' A more 
useful check comes from the positions of several radio sources 
found in Arp 299 ( |Neff et al.|20Q4] >. Two of those sources (Bl and 
D) match with X-ray sources and we find agreement to within half 
a pixel, with a hint of a small (0"25) offset of the X-ray image to 
the North. 

We also checked the astrometry of the optical discovery image 
(Newton & Puckett 2010) and the image taken by Joseph Brima- 
comb^Jand kindly provided to us. Comparison of the positions of 
15 USNO B1.0 and 2MASS stars gives an RMS scatter of 0"25. 
Together with the uncertainties in the USNO B1.0 and 2MASS po- 
sitions and the uncertainty of about 0"4 in the position of the super- 
nova on top of the galaxy, we estimate the accuracy of the optical 
position at 0"5. The best fit position of the supernova in these im- 
ages is RA = 11:28:33.811, DEC=+58:33:51.75, 0"4 away from 
the reported discovery position. 

Bon d et al.|p010| ) report observations with the WIYN 0.9m 
telescope that have been registered to the SDSS image, which then 
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Figure 2. Chandra images of the region around 2010O before the explosion. Top: observation 1641 (left) and observation 6227 (right). Bottom: same images 
Gaussian smoothed by 3 pixels. Our best optical position and the reported optical position of 2010O are indicated as the right two 0"5 radius circles. The 
SExtractor position of the source in the difference image (Fig.[3j is the left most circle. 



was registered to the archival Hubble Space Telescope ACS images. 
They report that the position of the supernova is consistent with 
a bright blue, slightly resolved object, likely a young star cluster. 
Fig. [T] shows the region of the supernova in the ACS image, with 
the different reported positions (including that of the X-ray source 
discussed below) indicated. The optical positions agree to within 
the uncertainties of 0"5. We have registered the ACS image with 
the single radio source D from |Neff et aL] ( |2004| ), leaving some 
uncertainty in the absolute calibration of the circles relative to the 
image. 



4 RESULTS OF THE CHANDRA ANALYSIS 

In Fig. [2] we show the two Chandra images of the region. WAVDE- 
TECT detects a source with formal significance of 9.2 sigma in ob- 
servation 6227 at position RA=11:28:33.813, DEC=+58:33:51.95, 
i.e. 0"15 from our best optical position and 0'.'5 from the reported 
discovery position. However, the X-ray source is extended and in- 
cludes a ridge of emission that is also present in observation 1641 
(see Fig.[2|. There is no detection of a source in observation 1641. 
However, at the position of the supernova there is clear diffuse 
emission (see Fig.[2|, complicating the analysis and artificially in- 
creasing the significance of the detection in observation 6227. We 
therefore subtracted image 6227 from (a scaled version of, see be- 
low) image 1641 and show the result in Fig. [3] There is a point 
like source that suggest there is a transient X-ray source at a posi- 
tion consistent with the position of the supernova. After the vari- 
able central source it is the most significant source in the differ- 
ence image. We ran SExtractor ( [Bertin &"~A rnouts 1996) on the 
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Figure 3. Colour image of the difference between Chandra observation 
6227 and the scaled 1641 observation. In red the soft (0.3-1 keV) in green 
the medium (1-2 keV) and in blue the hard (2-8 keV) X-ray image. The 
positions of sources found with WAVDETECT are indicated with yellow 
symbols for sources in observation 1641 and with white symbols for sources 
in observation 6227. The four sources found with SExtractor in the differ- 
ence image (see text) are indicated with red symbols. The reported position 
of 2010O and our best optical position are indicated as the two white 0"5 
radius circles. 

difference image and clearly detect the transient source at position 
RA=1 1:28:33.89, DEC=+58:33:51.6. 

In the following we proceed on the assumption that the local 
background has not changed significantly, i.e. that we can use ob- 
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Table 2. Relative sensitivity and counts in the diffuse emission around the 
position of the transient source in the two observations, expressed as ratio 
of observation 1641 to 6227. 



Band 


high ab. BB 


high ab. PL 


low ab. PL 


counts 


Low (0.3-1 keV) 


1.07 


0.91 


0.83 


1.18 


Middle (1-2 keV) 


1.61 


1.61 


1.61 


1.63 


High (2-8 keV) 


1.82 


1.80 


1.80 


1.70 


Full 


1.69 


1.42 


1.30 


1.45 



servation 1641 as an estimate of the background. One should keep 
in mind that it is quite possible that the source is also present in 
the earlier observation. Because the two observations were taken 
with different chips, we first determine the relative sensitivity to the 
diffuse emission for the two observations in three different bands 
(0.3-1, 1-2 and 2-8 keV) assuming: 

a) an absorbed black body of temperature 1 keV, taking the absorp- 
tion as Nh = 10 21 cm -2 , based on the reported reddening of the 
supernova spectrum of Ay = 0.65 ( [Mattila et al.|2Q10j >. 

b) an absorbed power-law with index 2 and 

c) a power-law with index 2 and only the expected foreground ab- 
sorption (Nh = 10 20 cm -2 ). 

We compare these to the detected counts in the two observa- 
tions in a region around the transient X-ray source with radius 100 
pixels. The results are shown in table[2] We conclude that the ratios 
measured in the two observations agree roughly with the expecta- 
tions (in particular for the absorbed black body). We thus divide the 
counts in observation 1641 by these ratios and assume that gives a 
reasonable estimate of the diffuse background. Within a 1"0 radius 
(which should contain 90 per cent of the flux for a source around 
0.4 arc min off-axis) we find 32 photons in observation 6227, while 
the background estimate based on the scaled observation 1641, is 
10.0. For standard Poisson statistics the likelihood of detecting 32 
photons or more for an expectation of 10.0 is 7 x 10 -9 , corre- 
sponding to more than 5.5 sigma. However, the background value 
is determined from only a small region in observation 1641, using 
14 counts, and thus has a relatively large intrinsic uncertainty. We 
therefore used a small Monte Carlo simulation, where we generated 
25 million backgrounds, taking the (Poisson) uncertainty on the 14 
counts in the background determination into account. The result 
is that the probability of finding 32 photons or more increases to 
3 x 10 -5 or about 4 sigma. 

Translating the counts in the difference image to X-ray flux 
using the black body model, we obtain an unabsorbed flux of 9.1 
xlO -15 erg cm -2 s -1 . 

The counts in the different bands (0.3-1, 1-2 and 2-8 keV) for 
the difference between the scaled observation 1641 and 6227 are 
8.9, 11.3 and 1.8 respectively. This is softer than the spectra used 
above to determine the overall sensitivity difference. Indeed, the 
counts are consistent with a black body of temperature of 0.225 
keV 

We thus repeated the calculation of the fluxes with a 0.225 
keV black body and find an absorbed flux of 8.6 x 10 -15 erg cm -2 
s -1 and an unabsorbed flux of 1.2 xlO -14 erg cm -2 s -1 . For a 
distance of 41 Mpc this corresponds to an unabsorbed luminos- 
ity of 2.5 xlO 39 erg s -1 . A Monte Carlo calculation generating 
the background as described above and adding sources of different 
luminosity provides the 1 sigma uncertainties on this luminosity: 
L x = 2.5± 0.7 x 10 39 ergs -1 . 



5 DISCUSSION AND CONCLUSION 

If the X-ray transient found in the second Chandra image is asso- 
ciated with the SN2010O (see below for a discussion of the possi- 
bility of a chance alignment), it suggests that this is the explosion 
of the second star in a binary system. Before the supernova explo- 
sion the system must have consisted of a Wolf-Rayet star orbiting a 
compact object. We can compare the X-ray properties of the tran- 
sients to the outbursts seen in Galactic black-hole X-ray binaries 
(see McClintock & Remillard 2005 for an overview). The peak lu- 
minosity is higher than those of the Galactic sources, but these are 
particularly poorly known, as the distances to these systems are not 
well determined (see Jonker & Nelemans 2004 ). The outburst lu- 
minosity is below the Eddington limit of a 10 Mq black hole, if we 
assume pure helium accretion. The spectrum is much softer than the 
typical 1 keV found in Galactic black hole binaries. However, the 
very soft spectrum and high luminosity fit very well in between the 
highest luminosity point of the Galactic X-ray binary XTE J1550- 
564 and the ultra-luminous X-ray binaries NGC 1313 X-l and X-2 
(see figure 7 of |Soria|2007) . 

The transient thus could be a relative of the nearby Wolf- 
Rayet X-ray binaries Cyg X -3 (van Kerkwijk et al71fT996] >, IC 10 
X-l jClark & Crowther|2004||Bauer & Brandt|2004l > and NGC 300 
X-l (Carpano et al. 2007) and the soft spectrum suggests a (mas- 
sive) black hole accretor. However, all of these seem to be fairly 
persistent X-ray sources, with luminosities well below the Edding- 
ton limit. The transient thus would represent the equivalent of the 
outburst in soft X-ray transients, in a Wolf-Rayet X-ray binary. 

With a distance modulus of 33 and and a reported visual mag- 
nitude of 15.8 around the peak (Newton & Puckett 2010), together 
with the extinction of 0.65 as derived from the spectrum (Mattila 
|et al.|2010) , the absolute magnitude of SN2010O is -17.9. This is 
well within the large range of absolute magnitudes of type lb super- 
novae (Richard son et al.|2006| ). Unfortunately, there are no known 
correlations of lb supernova explosion features and the properties 
of the progenitor Wolf-Rayet stars. 

The whole discussion above depends on the association of the 
X-ray transient with the supernova. With WAVDETECT 22 unique 
sources are found within half an arc min from the position of the 
supernova. If we assume we would consider sources within \" as a 
possible match that would give a chance alignment probability of 2 
per cent. Ten of these sources are actually within 15" and the differ- 
ence image in Fig. [3] shows two sources not detected by WAVDE- 
TECT if we include the transient. We therefore conclude that there 
is a 5 per cent probability of a chance alignment. However, the close 
match in position and general consistency of the X-ray properties 
with what could be expected from an X-ray binary in which the 
second star exploded as a lb supernova, make the connection plau- 
sible. The possible association of the supernova with the cluster 
seen in the HST ACS observations (Bond et al. 2010) also allows 
for the possibility that the X-ray source is coincident with the po- 
sition of the SN, but is actually an unrelated X-ray binary in the 
cluster. However this is not very likely, as there are only 1 8 X-ray 
sources found in Arp 299, and many more star clusters, and a study 
by |Kaaret et al.| ( [2004| shows that in starburst galaxies the clusters 
have only rarely an X-ray source right on top of them. 

We therefore conclude that the observations presented in this 
paper suggest the exciting possibility that supernova 201 0O was 
the explosion of a massive Wolf-Rayet star that was part of a X-ray 
binary containing a black hole. 

Further study, in particular for other supernovae with pre- 
supernova Chandra data, will be needed to address the question 
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of how often type lb and Ic supernovae explode in X-ray binaries. 
We are in the process of analysing the earlier type lb and Ic super- 
novae for which there is archival Chandra data (Nielsen et al. in 
preparation). 
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